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1. Introduction 

The decoherence process followed by a given physical system is usually modelled by a 
(generalized) master equation for its density operator. In the literature of open quantum 
systems, master equations are usually called Markovian if they did not involve an explicit 
time integration. Otherwise, they are termed non-Markovian. These designations are 
employed in this paper, albeit they do not agree with the mathematical definitions 
of Markovian and non-Markovian processes, as was shown long ago pQ. Markovian 
master equations, which can be derived from system-environment models employing 
the Born-Markov approximation [2], sometimes present unphysical behavior, because 
they lead to negative density operators [31 0]. Master equation of the Lindblad form 
[5j [6] , which preserve the fundamental properties of density operators (non- negativity, 
unit trace and hermicity), are also markovian. Under special conditions they can also 
exhibit unphysical behavior, like the unbounded increase of the energy of the system 

0- 

Markovian master equations are important tools to model open quantum systems, 
not only due to their mathematical simplicity, but also because they capture the physical 
behavior of many important systems, such as open QED (Quantum Electrodynamics) 
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systems [El El QUI EEE]- Many Markovian master equations are obtained using the Born- 
Markov approximation, which relies on a series of assumptions which are not always 
satisfied p2] : weak system-environment coupling, separable total density operator, bath 
correlation time much smaller than the relaxation time of the system, and unperturbed 
transition times of the system much smaller than its relaxation time. The next to last 
assumption implies that the environment is infinite. Thus, loosely speaking, Markovian 
behavior is typical of systems which interact weakly with infinite environments. 

One of the alternatives to the use of generalized master equations to describe 
open quantum systems is the Feynman- Vernon influence functional [13]. Other, which 
was derived by Kraus [2] employing the idea of complete positivity, is known as the 
"operator-sum representation" of the quantum dynamics of an open system. The 
operator sum representation and Lindblad-form master equations can be related. In 
effect, if a closed quantum system S, initially in a product state ps(0) = pa(0) <g> ps(0), 
comprises two interacting subsystems, the system of interest A and its environment B, 
the exact dynamics of the state any of the subsystems, can be put in Kraus form. If the 
inequality r c <C r <C r H is satisfied, where l/r c is the cutoff frequency of the bath density 
of states, t an adequate coarse-graining time scale and th the characteristic time-scale 
of the hamiltonian evolution of the system, then it is possible to derive a completely- 
positive master equation starting from the Kraus operator-sum representation JT5J [16] , 
linking both descriptions. The operator-sum representation, used to represent noisy 
dynamics, is popular in the area of quantum information processing. In particular, 
Nielsen [TF] introduced several archetypal quantum operations on a qubit, including the 
depolarizing, phase-damping (or phase-flip) and amplitude-damping channels, which are 
customarily assumed Markovian in theoretical and experimental analyses [18] . However, 
the description of decoherence in solid state systems, the most promising scalable 
realizations of quantum processors, often need to be non-Markovian [191 EH [ZD 122] - 

The purpose of this manuscript is to show, explicitly, that the archetypal one-qubit 
noisy channels generally describe non-Markovian dynamics. After a brief review of the 
three channels that we consider in this paper, we present simple (yet non-Markovian) 
microscopic models for them. The decoherence process is caused by the fluctuation 
of macroscopic variables which enter in the Hamiltonian or by the establishment of 
correlation betwen the system and its environment. 

The standard form of Lindblad master equations [5] for the density operator of the 
system of interest, p, is 

dp(t) i 



dt h 



where H is a Hamiltonian and {Lj} are (possibly non-hermitian) operators, usually 
known as Lindblad operators. For a given physical process the decomposition (Tj[|) is not 
unique. However, it is convenient to interpret the first term of the right hand side of 
(TjQ) as the unitary evolution in the absence of interaction with the environment, and 
the second term as the environment-system coupling. On the other hand, the Kraus 
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representation 

p(t) = ^^(t)p(0)4 t (t), (2) 



where the Kraus operators Ei(t) satisfy the condition of probability conservation 
E\(t) Ei(t) = I, gives the evolution of the density matrix of the system of interest. 
We consider a qubit with orthonormal states |0) and |1). In this basis the 
operator a 3 = |0) (0| — |1) (1| is diagonal, and the lowering operator is given by 
<t_ = |1) (0| = Si — i<72 = o-+. Several noisy channels have been defined for a single 
qubit. For example, the depolarizing channel 

/3(p)=P^ + (l-p)p(0), (3) 

describes a process where the state of the qubit, initially in its state p(0), remains in this 
state with probability p, and changes to the maximally mixed state 7/2 with probability 
p. If a two-qubit mixed state is used, instead of a Bell state, in the standard teleportation 
protocol, a generalized depolarizing channel is obtained [23]. An amplitude damping 
channel is obtained when a qubit interacts with a large reservoir at zero temperature, 
or a two-level atom in the electromagnetic vacuum emits spontaneously [21 |8l 121 Ell ED] ■ 
An exciton confined to a quantum dot and coupled to phonons [24J, like a spin in 
a magnetic field whose strength fluctuates in time, suffers a dephasing process. The 
Markovian master equations corresponding to the depolarizing, amplitude-damping and 
phase-damping channels are 



i 

(2 <7_ p(t) a+ - + a_ p(t) - p(t) a + 0_) , (5) 



alt 

dp(t) _ r 

dt ~ 2 
dp(t) 



T (2 p(t) a + a_ - &+&-p(t) - p(t) <r + a_) , (6) 

respectively. The operator-sum representation (T5]) of these channels depends on the 
Kraus operators 



E o=\l 1 --^ 1 ' E i = ^-a i , i = 1,2,3; (7) 



4 ' 2 

| =|0)(0| + 7T^|1)(1| andE 1 =^Vp|0)(l|; (8) 
E = ^/T^I, Ei = ^p\0){0\, E 2 = y/p\l){l\- (9) 

where p(t) = 1 — e _rt . A given physical process can be described by multiple Kraus 
representations, as can be seen in the phase-damping channel ([9j), which can also be 
given by the two Kraus operators [25J 



£ = Ji-f/, £i = Jf*s. (10) 

In the following sections we show that the same Kraus representation can be related to 
many master equations, generally non-Markovian. 
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2. Simple quantum microscopic model of a depolarizing channel 

Electronic spins in semiconductor quantum dots have been proposed as a prototype of a 
scalable quantum computer [26]. A Caldeira-Leggett microscopic model of decoherence, 
considered to be a generic phenomenological description of the environment of the spins, 
have also been advanced. The interaction of the spin with the collection of harmonic 
oscillators, which describe its environment, is of the form 

Him = A ^ a i 9k ( aik + ^k) ' 

i k 

where (aj fc ) is the annihilation (creation) operator of the k — th harmonic oscillator 
coupled to the i—th component of the spin. The interaction is clearly isotropic. It is also 
natural to represent the degrees of freedom of the environment as spins, corresponding 
to the hyperfme interaction, or to the interaction of the electronic spins with other 
electronic spins. Here we consider a simple variant of the spin star model [271 [28] in 
which the dynamics of the spin-bath system is governed by the total Hamiltonian 

N 
k=l 

where = ho~i/2 is the i — th component of the electronic spin and Ski is the i — th 
component of the k — th environmental degree of freedom, which, for the sake of 
simplicity is assumed to be a one-half spin. The coupling constant g has units of 
frequency. It is convenient to write the environmental states in the basis of collective 
angular momentum. It is well-known that the bath states can be labeled as \l,mi) 
where n < I < N/2, mi = —I, — I + 1, • • ■ , I and k = (k = 1/2) if N is even (odd). The 
degeneracy of the bath momentum I is given by J27J [29], [28] 

"W> = (^ / )-(f_ J L 1 )- (13) 

In the new basis the spin-bath Hamiltonian is simplified to H — ®zfs ■ Si, where the 
degeneracy must be taken into account. In the variant we consider here, described by 
the Hamiltonian 

H = ®i^s-S h (14) 

h 

not only the degeneracy is not given by (|T3|) but the coupling coefficients can be different. 

The initial state of the total system is assumed to be factorized pr(0) = 
As(0) <S> P_b(0), that is, equal to the product of the initial spin state ps(0) and the 
initial bath state, ps(0). The latter is assumed to be the maximally mixed state 
2~ n In x n. The notation I^xn stands for the identity in N dimensions. In the new 
basis we have the interaction between spin 1/2 and spin /, where the initial state is 
ps(0) ® {21 + l) _1 I(2i+i)x(2/+i)- For a fixed value of I it is convenient to make a unitary 
transformation, from the separated basis, whose states ||, /; m s , nij) are simultaneous 
eigenstates of s 2 , Sf, s z and J z , to the coupled basis, whose states \\,l;j,mj) are 
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simultaneous eigenstates of s 2 , Sf, J 2 and J z , where J = s + S[. In terms of these 
operators the spin-bath Hamiltonian reads ©/|(J 2 — s 2 — Sf). One can go from the 



separated basis to the coupled one using 



2, I] I — 2, I + 2 — P) 



cos ov 



sin a 



— sin a p cos a p 




i,Z;-i,Z-p+l> 
| 2, Z; 2? Z — p) 



(15) 



where the elements of the transition matrix are given by the Clebsh-Gordan coefficients 



cos ct r 



V 



sma r 



2l-p+l 



"~ + V p ~ V 21 + 1 

The dynamical problem is solved, for the initial condition 



lV>o) 



1 1 



'2' 2 



1, 



) (g) |Z,mj) . 



1 

2' 2' 

If the initial state is expanded in the coupled basis, evolve the state taking into account 
that the eigenenergies hgl/A and — hg(l + l)/4 correspond to the eigenstates of Sf with 
eigenvalues / + \ and Z — |, respectively, and transformed back to the separated basis 

2' ' 2 ! 

i Z ._I 

2' ' 2 ! 



mi 



-^,77li + 1 > + C™' 



1 / 1 



+ c 



mi 



1 1 

2' Z; 2' m ' _1 



(16) 



The time dependent coefficients c" 



mi 



1, 2, 3, 4 are given by 



-2ic+e 4 



v /(/-m / )(/ + m z + l) . f(2l + l)gt 



c 2 



c 4 



c + e l * 



— 2ic_e 



cos 



21 + 1 
(2/ + l)gt 



sin 



2m; 



2/ + 1 



cos 



(2Z + 1)^ + ^2m. 



2/ + 1 



sin 



sin 



(21 + ljgf 
4 

(2/ + 1)^ 



V^T+miKZ -mi + 1) . f(2l + l)gt 



(17) 
(18) 
(19) 

2/ + 1 -i ' ' < 20 > 

After finding the total density matrix operator, and tracing out the bath degrees of 
freedom, the reduced density operator for the spin one-half system is found to be 

p mi = [\c?f + \cT\ 2 ) |o> (o| + cT^\o) (i| + c7^\\) (o| + (Kf + \cT?) |i> (i| , 

where the state ||, — |) (||, |)) have been identified with the state |0) (|1)). Now, it is 
possible to take into account that the initial bath spin-/ state is (2j + l) _1 I(2j+i)x(2j+i) ) 
to show that 



sin 



Pi 



21 + 



1 1 

— y Pm = -(hx2 + Si(t) 



(21) 



can be calculated exactly. The Bloch vector si(t) completely characterizes the quantum 
state of spin one-half system and is given by 

4Z 2 + 4/ + 3 + 81(1 + 1) cos ( ^M^j 



si(t) = 8i(oy 



3(21 + l) 2 



(22) 
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where Sj(O) is the value of the Bloch vector of the initial state of the system. For large 
values of I, si(t) ~ sj(0)(l + 2 cos ( ^• 2l+ ^ l9t j)/3. If we have a fixed value of I and a 
gaussian distribution of coupling constants, in which the probability density to have a 
particular value g for the coupling constant is 

p G (g) = (27ra 2 )- 1/2 exp(-//(2a 2 )), (23) 
a being the standard deviation of the distribution, then the average Bloch vector 

*(«) = «(0) 3^ , (24) 

decreases until its length is reduced for a factor x = f^^nr , which satisfies the 
inequality | < x < |. 

If the depolarizing channel ([7]) is applied to an initial pure state, its Bloch's vector 
evolves as so(t) = Sd(0)(1 — p). We thereby conclude that we have a depolarizing 
channel for which 

p(i) = |lTTF (1 "^ <2 ' +1,W/8) ' (25) 

a gaussian decay very different from the Markovian process, p(t) = 1 — e~ rt . Moreover, 
in constrast with the Markovian case, where p(t) varies from zero to one, in the 
non-Markovian process described by ( l24l) . p(t) varies, monotonically, from zero to 
Pinf(0 = 3(2i+i)2 < §, i-e., this process does not display complete depolarization. 

As a second example of the distribution of coupling constants, we assume a 
Lorentzian distribution p(g) = a(Tc(g 2 + a 2 )) -1 , where a is the scale parameter which 
specifies the half-width at half-maximum. The average Bloch vector decreases as in the 
previous example 

, N 4Z 2 + 4/ + 3 + 81(1 + l) e -( 2i+1 ) a */ 2 , . 

Sl (t) = sM , (26) 

but in an exponential way. The time-dependence of the parameter p, given by 

P(t) = 3^W (1 - e_(2/+lW/2) ' (27) 
is similar, but different, to that of a Markovian process. Indeed, the corresponding 
master equation reads 

= ^ S (2 ^ &i ~ &i &i ~ p® bi ^ ' (28) 

which does not have the Lindblad form because now ^(t) is an explicit function of time 
_ / 4/ 2 + 4/ + 3 \ 

7(*) - r \i - — 2 - 4/ - 3 - g/(/ - 1)e _ (2m)at/2 J ■ (29) 

In a simpler, but also interesting case, / and g fixed, p(t) varies periodically, and 
the master equation is given by (1251) with a time- dependent 7(t), 

/ (2l+l)gt\ 



(21 + l)g 8/ ( / + !) sin 



t(*) = v 9 ' ;, 2/+1wV (30) 

2 4Z 2 + 4Z + 3 + 8Z(/ + l)cos(^±^J 
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which, in contrast to (1291) . attain negative values. 



s(t) = s(0) oos(2^) + £ x S (0)^^ + (* ■ s(0)) ^ , 2 & ^ , (33) 



3. Simple classical microscopic model of a depolarizing channel 

In nuclear magnetic resonance experiments and Bose-Einstein condensates the 
decoherence process is often caused by residual fluctuating magnetic fields, see for 
example [30], which can be considered classical. In the Hamiltonian description of 
this process 

H = -n- B = hg£- <r, (31) 

the dimensionless independent random variables i — 1, 2, 3, assumed to be gaussian of 
zero average and standard deviation a, are proportional to the corresponding magnetic 
field components. The constant g have units of frequency. If the initial state of the spin 
system is p(0), its state at time t is 

p{t) = >(0)e^ CT = ~ (I 2x2 + s(t) ■ a) , (32) 

where 

, sin(2^0 + ( £ . „ (0)) $(l-cos (2fftQ) 

and £ = V£ ' Averaging over the different realizations of the noise variables £j, 
z = l,2, 3, the state of the two-level system is 

m = \ ( J 2x2 + «(o • v) , 8i {t) = 8i {o) (™^W) + ^(i-cos(^e)) j = s . (0)/(i)) (34) 

where the propety of vanishing averages of the distributions of £j were used. Employing 
the gaussian probability distribution pd^i), i = 1,2,3, fl23l) . and transforming to 
spherical coordinates given by r = £, 9 = arccos(£s/£), = arctan(£ 2 /£i), where 
r G [0, oo), G [0, 27r) and 6 G [0, tt], one finds for the polarization factor f(t) 

/(*) = \(l + 2(1 -AgW)e^ W ). (35) 

In figure [1] we have plotted the polarization factor using o = 1 and time measured in 
units of 1/g. 

We stress that albeit none of the examples considered in the previous section and 
this section are Markovian, all of them are described by the Kraus representation (J7J) 
of the depolarization channel, with different definitions of the parameter p. The time- 
dependence of p allows for non-exponential decay and even recoherence, in which case 
7(t) is negative. 

4. Simple quantum microscopic model of a dephasing channel 

The model hamiltonian we consider, given by 

frwjt) v^t ( , , Cka A fir\ 
H = -^—a z + > ftw k [a k ^ a fe H , (36) 
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Figure 1. Polarization factor f(t) for a microscopic model with classical gaussian 
noise, with zero average and unit standard deviation. The continuous line represents 
the theory and the dashed line a numerical average over 10000 realizations. The time 
t is measured in units of 1/g. 



describes the interaction of a two-level system with a collection of oscillators, through an 
interaction which conserves the system's observable a z . The creation and annihilation 
operators a\ and a k satisfy the usual boson commutation relations [a^aj.,] = 8 kk >. 
Hamiltonians of the type we consider here have been explored by many authors 
pTTj [321 E31 EH ES]- We assume an initial state of the Feynman- Vernon form p(0) = 
Ps(0) <8> rifcPfe; where p k = exp ^— j3hw k a\a k j /Z k is the thermal state of the k — th 

mode at inverse temperature (3 = l/{k B T) and Z k =tr ^exp ^— f3hu k ala k )^ . If the 
action of the unperturbed Hamiltonians is separated from the interaction, the total 
density operator is 

p{t) = e-^ Waz/2 e- iE ^ fc ^ a W / (t,0)p(0)W / t (t,0)e iE ^' i ^ a e-^ ( ' ),Tz/2 (37) 
where the interaction evolution operator in the interaction picture, Ui(t,0), 



Ui(t,0) = Texp l-i / c/rV (c k e ibJkT a\o z + h.c.) = Te~ l ] fo d ^ k H k (r) ^ 




needs the time-ordering prescription, indicated by r, because the commutators 
[H k (r), H k (r')} do not vanish. However, due to the simple Lie algebra satisfied by 
the operators appearing in the interaction Hamiltonian, [a k o~ z , a\a z ] = 8 kl i, and tracing 
out the environmental degrees of freedom we can write 

ps (t) = e- in{t) ^ /2 Ti k ^e-^o^E fc ^(r) ps(0 ) JJ pke -i!ZdTZ k H k (T)j e -m)^l\ (39) 

Employing algebraic techniques [2] to calculate the trace, the density operator of the 
system can be written as 

p s (t) = e ~^W 2 ( e - 2r ( t )( 1 -^^)p s (0)) e' in ^/ 2 , (40) 

where the functions Q(t) = J Q *d™(r) andT(t) = J2 k \c k \ 2 (l—cos(u k t)) coth(hui k f3 /2)/cul 
have been defined. Taking into account that ps(0) = Ylij Pij N) 01 the expression into 
brackets in equation (j4"0"l) can be simplified as follows 

e-arWd— -) P5 ( ) = e -2r(*)(i— . ) £ p .. {i) {j{ = £p.. e - 2 r(*)(i-^>*(;)) (i> (j | , (41) 
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with s 2 (0) = 1 = — Sje(l). If we explicitly write the four terms we have 
c -ar<t)<i-*.-T.) p5 ( ) = poQ | ) ( | + e- 4r Wp i |0> (1| + e^Vo |1> (0| + Pu |1> (1| , 
which can be recast as 

e - 2 r (i )(i-. z . CTz)ps(0) = e -4r Wp(Q) + (1 _ e - 4 r( t)) (poQ |Q) (Q| + pn ^ (42) 

We recognize the form of a phase damping channel with p{t) = 1 — e~ r ^\ Finally, 
including the effect of the unitary operators e ±lQ ( t )°'^/ 2 ^ we se e that, at time t, the 
reduced density of the system can be written as a phase damping (PD) channel, with 
Ei = e- iQaz/2 Ei, and E = y/T=pI, E 1 = y/p |0) (0|, and E 2 = y/p |1) (1|, as before. 

The behavior of p(t) depends on the number of oscillators of the environment, their 
frequencies, coupling constants and the temperature of the bath. If the "environment" 
contains a single oscillator, p varies periodically, with the unperturbed frequency of this 
oscillator, between and pu'- the higher the temperature (or the stronger the coupling 
or the smaller the frequency) the greater the value of pu- In this case, p{t) does not 
have a limit for long times. Moreover the master equation, which can be written as 

^ = -^[hw(t)a z , p] + 7 (t) (2a z pa z - a 2 z p - pa 2 z ) (43) 

is non-Markovian. In particular, j(t) is negative for an infinite number of time intervals. 
When the limit to the continuum is taken (J2k ~^ I doj/(2ii)), in the decoherence 
function T(t) at finite temperature a zero-temperature contribution can be isolated, 

r (t)= r^^(l-coB(ut))du. (44) 

The behavior of the remaining contribution Tp(t) = T(t) — To(t), often differs 
significantly at large times from that of r (t), as shown in Figure [2J The examples 
given before illustrate how p{t) can oscillate, or grow to one (or to a constant smaller 
than one) exponentially or not. 




01 2345678 
t/T 



Figure 2. Dephasing channel. Function p{t) for a single oscillator with frequency 
1/(2T) and |c| 2 coth(fiw/3/2)/w 2 ) = 2 (solid line), and for a continuum of oscillators 
with D(ui) — 27rwe _LjT for zero temperature (dashed line) and f3 = 1/T (dotted line). 
Time is measured in units of T. 
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5. Simple classical microscopic model of a dephasing channel 

A dephasing channel also occurs when the magnitude of a classical magnetic field 
changeas randomly. The Hamiltonian which describes this situation, H = hw(t)a z /2 + 
hg£(t)a z , where is a random variable. The equation of motion for the density 
operator is easily integrated, 

p(t) = e -m)i°*, ■] e -^Lfo^V--] p (0) 

= e - im \°*A tist f dtl . . . fduik) ■ ■■at n )[a z , -]>(o), 

where fl(t) = ^u{t)oIt. A gaussian stationary process with zero average satisfies 
[36] 



e(ti)e(t 2 ) =$(t 1 -t 2 ), 

f(*l)"-f(*2n-l) =0, 

■ ■ ■ £(t 2n ) = ~~ *«) ' ' ' $ (^(2n-l) - *i(2n)), 

Perm 

where n is apositive integer and the overline indicates the expected value. The averaged 
density operator 

00 < ;„\in rt rt 



( ' \2n pt {•% 



n=0 

X 

Perm 

can be simplified taking into account that there are (2n)\/(2 n n\) permutations of the 2n 
times ti 

m - t Qf*jfW - «.))" [*.. ■]>(»)■ 

The dynamics of the two-level system 

p/ t \ _ e -ai(*)k..-] e -iJ?*iJ?«tta*(*i-*a)k.,k..-]]p(o), (45) 

corresponds to a phase damping channel with = 1 — e -1 ^, where T(t) = 

\ jldti j^dt 2 ^{t\ — t 2 ), which generally is non-Markovian. Only if the two-point 
correlation function $(ti — t 2 ) = o- 2 5(t\ — t 2 ) the dynamics is Markovian. 



6. Simple quantum microscopic model of an amplitude damping channel 

Finally, we consider a qubit interacting with a collection of harmonic oscillators, 
which models the degrees of freedom of its environment, described by the Hamiltonian 
H = r f |1) (1| + J2 k Mfc4 a * + Efc {ckalo-- + 4a fc <x+) [33 H21 EH]. Here, a k is the 
annihilation operator of the k — th mode of the environment of free frequency u k , and c k 
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are interaction constants. If the initial state of qubit and environment is \ip) s <8> Ylk |0)fc> 
at time t it can be written as 

i*(t)> = ait) ii) ® n io) fc + p(t) io) ® n io) fe + i°) ® ® n i°>* • ( 46 ) 

ft k l k^l 

The use of Schrddinger equation leads to a set of coupled first-order linear equations for 
the coefficients a(t),(3(t) and ji(t). Solving the equations for "fi(t) in terms of a(t), and 
replacing into the equation for a we find the integro-differential equation 

^P- + iua{t)+ f dry2Wi\ 2 e- iu>l{t - T) a(r) = 0. (47) 
dt J 

The solution of this equation, which is also obtained for the independent oscillator model 
[39] . can be written as a(t) = a(0) exp(— A(t) — iQ(t)), where A(i) > and A(i) are 
real functions. The reduced density of the qubit is obtained tracing out the state of the 
oscillators 

p s (t) = tr E |*(t)> = fi(t) |1) (1| + (1-M*)) |0> (0| + (a(t)/3*(0) |1) (0| + /i.e.) , (48) 

where /i(t) = |a(£)| 2 = |a(0)| 2 exp(— 2A(£)). A unitary contribution to the dynamics of 
the qubit can be isolated, 

p s (t) = e~ iU ^ (/i |1) (1| + (1 - /i) |0) (0| + (a(t)/3*(0) |1) (0| + h.c.)) e^ 1 ^. (49) 

Now, taking into account that 

|0)(0| + V^pllXll) (a |1> + P |0» (a* (1| + /T (0|) (|0)(0| + v^^UI 

= (l-p)|a| 2 |l)(l| + |/3| 2 |0)(0| + (v/l-pa/3*|l) (0| + /i.e.) , and (50) 

y/p |0)(1| (a|l> + /? |0» (a* (1| + /T (0|) ^ |1>(0| = p|a| 2 |0)(0| , (51) 

where |a| 2 + \/3\ 2 = 1, we conclude that p(t) = 1 — e~ 2A( ^. Moreover, we see that the 
qubit dynamics corresponds to an amplitude amplitude damping (AD) channel with 
Kraus operators E, t = e'^^Ei, and E = |0}(0| + y/T=p |1)(1| and E x = y/p\0)(l\. 
This amplitude damping channel is, as the noisy channels considered in the previous 
sections, generally non-Markovian and includes cases in which the environment consists 
of a finite oscillators (even one). 
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